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Abstract 

We propose a model addressing the coupling mechanism between two spatially modulated mono- 
layers. We obtain the mean- field phase diagrams of coupled bilayers when the two monolayers have 
the same preferred modulation wavelength. Various combinations of the monolayer modulated 
phases are obtained and their relative stability is calculated. Due to the coupling, a spatial modu- 
lation in one of the monolayers induces a similar periodic structure in the second one. We have also 
performed numerical simulations for the case when the two monolayers have different modulation 
wavelengths. Complex patterns may arise from the frustration between the two incommensurate 
but annealed structures. 
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I. INTRODUCTION 



Quite a number of physical, chemical, and biological systems manifest some type of 
modulation in their spatial ordering Such structures are stripes and bubbles in 

two-dimensional (2D) systems, or lamellae, hexagonally packed cylinders, and cubic arrays 
of spheres in three-dimensional (3D) CcLSGS cLS well as more complex structures such as gy- 
roids. Examples of such systems include ferromagnetic layers [4], magnetic garnet films |5|, 



ferrofluids 
copolymers 
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£], dipolar Langmuir films ItI, rippled phases in lipid bilayers and block 



10|] . Modulated phases may also occur in systems described by two (or more) 



coupled order parameters, each favoring a different equilibrium state [llj. The observed 
spatial patterns exhibit striking similarity even for systems that are very different in their 
nature. It is generally understood that the modulated structures are formed spontaneously 
due to the competition between short- and long-range interactions. 

In the case of 2D ferromagnetic layers, for example, the short-range interaction arises from 
magnetic domain wall energy, while the long-range interaction is due to magnetic dipole- 
dipole interaction which induces a demagnetizing field [2], sl. Adding both contributions 
and minimizing the total free energy with respect to the wavenumber g, one obtains the 
most stable mode q* 7^ 0. This description is valid in the weak segregation limit (close to a 
critical point), where the equilibrium domain size is given by d* = 27r/g*. In general, this 
quantity depends on temperature and/or other external fields. 

In this paper, we shall consider two modulated monolayers that are jointly coupled. Our 



motivation is related to recent experiments 



3y Collins and Keller 1^ who investigated 



Montal-Mueller planar bilayer membranes [13|] composed of lipids and cholesterol. In this 
technique, a bilayer is constructed by preparing separately two independent monolayers and 
then combining them into one joint bilayer across a hole at the air/water interface. The 
experiments addressed specifically the question of liquid domains in the two leaflets, and 
the mutual influence of the monolayers in terms of their domain phases. In the experi- 
ment, asymmetric bilayers are prepared in such a way that one leaflet's composition would 
phase-separate in a symmetric bilayer and the other's would not. In some cases, one leaflet 
may induce phase separation in the other leaflet, whereas in other cases, the second leaflet 
suppresses domain formation in the original leaflet. These results imply that the two leaflet 
coupling is important ingredient in determining the bilayer phase state. 
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Motivated by these experiments, the coupled bilayer system was investigated theoretically. 
The coupling mechanism arises through interactions between lipid tails across the bilayer 
midplane, and the phase behavior of such a bilay;er membrane was computed using either 
regular solution theory ij] or Landau theory 15|]. The theoretical results are in accord 
with several of the experimental observations. It should be noted that all previous models 
dealt with the coupling between two macro-phase separated leaflets, while it is also of 
interest to investigate the coupling between two micro-phase separated (modulated) leaflets. 
Furthermore, one might also consider the interplay between a macro- and a micro-phase 
separation. 

In the present work, we suggest a model describing the coupling between two modulated 
systems, and, in particular, we analyze the influence of this coupling on the phase behavior of 
two coupled 2D monolayers. When the two monolayers have the same preferred periodicity 
of modulation, we obtain the mean-field phase diagrams which exhibit various combinations 
of micro-phase separated structures. In some cases, the periodic structure in one of the 
monolayers will induce a modulation in the other monolayer. Interesting situations take 
place when the two monolayers have different preferred wavelengths of modulation. Here 
the frustrations between the two competing modulated structures need to be optimized. 
These structures and their dynamical behavior are examined using numerical simulations. 
Although there has been so far no experiment which directly corresponds to the proposed 
model, our predictions may be verified, for example, by constructing Montal-Mueller bilay- 



point 



ers 13 1 out of two lipid monolayers that exhibit a striped phase near the miscibility critical 
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17|. 



In the next section, we present a phenomenological model describing the coupling between 
two modulated lipid monolayers. In Sec. IIIIl we discuss the case when the two monolayers 
have the same preferred wavelength of modulation. Monolayers having different preferred 
wavelengths will be considered in Sec. HVj and some related situations are further discussed 
in Sec. |Vl Although we limit our present analysis to 2D systems, the suggested model can 
be generalized to 3D systems as well. 
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FIG. 1: Schematic illustration of two coupled modulated monolayers forming a bilayer membrane. 
Each monolayer is composed of a binary A/B lipid mixture, which can have a spatial modulation. 
The relative composition of the two lipids in the upper and the lower leaflets are defined by (p and ip, 
respectively. In general, the average composition in the two monolayers can be different. The lipid tails 
interact across the bilayer midplane. The phenomenological coupling term between these two variables 
are assumed to be bilinear of the form —Acptp in the free energy of Eq. ([T]). 

II. MODEL 



In order to illustrate the coupling effect between two modulated systems, we imagine 
a pair of lipid monolayers forming a coupled bilayer. Each of the monolayer can undergo 
separately a micro-phase separation. As shown in Fig. [H we assume that each monolayer is a 
mixture of two lipid species, say lipid A and lipid B. Their area fractions are defined by 0a (r) 
and 0B(r), where r = {x,y) is the 2D positional vector. By assuming that the monolayer 
is incompressible, 0a (r) + 0B(r) = 1, the monolayer composition can be characterized by a 
single order parameter defined by the relative A/B composition 0(r) = 0a (r) — 0b (r)- Let 
us denote this local order parameter of the upper and lower monolayers by 0(r) and tp{r), 
respectively. The coarse-grained free-energy functional for the coupled modulated bilayer is 
written as: 



F[0,V^] 



F„[0] + F, 
dr 



A / dr. 



2(VV)^ - 2(V0)^ + ^0^ + \<P' 



72„;.\2 o^/r7„;.N2 , ''"„/.2 , ^ -.4 



+ 2D{V'^y - 2C{V^y + -V + - /iv^ - A0^ 



(1) 



This is a modified Ginzburg- Landau free energy expanded in powers of the order parameters 
and -0 and their derivatives. The -Fu[0] free energy has five terms depending only on 
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and its derivatives. It describes the upper monolayer and its possible modulations, while 
the coefficients of the Laplacian squared, the gradient squared and the 0^ terms are taken 
to be numbers, for simplicity. Similarly, F([^] describing the lower monolayer contains the 
next five terms that are only functions of ip and its derivatives. The last term represents 
the coupling between the two leaflets as will be explained later. The coefficients of the two 
gradient squared terms are both negative (C > 0) favoring spatial modulations, whereas the 
coefficients of the Laplacian squared terms are positive {D > 0) to have a stable modulation 
at finite wavenumbers. The 0^, 0^, ip"^ and ■0^ terms in F are the usual Landau expansion 
terms with r = (T — Tc)/Tc being the reduced temperature (Tc is the critical temperature). 
For simplicity, the two leaflets are taken to have the same critical temperature (and hence 
the same r). Finally, the linear term coefficients, fi^f, and fi^, are the chemical potentials 
which regulate the average values of and ip, respectively. 

In the absence of the coupling term (A = 0), each of the two leaflets can have its own 
modulated phase. Free energy functionals such as Fu have been used successfully in the past 
to describe a varietyof modulated systems: magnetic garnet films p] , Langmuir films [7| , 
diblock copolymers jol. . and amphiphilic systems 19]. Furthermore, interfacial properties 
between different coexisting phases have been investigated using a similar model 
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mm 



In the above expression for the free energy F, the 0-leaflet has a dominant wavenumber 



= 1/ \/2, and so has the -^-leaflet with = ^JC jlD. The modulation wavenumbers and 
amplitudes of the two monolayers coincide when D = C = \ and the average compositions 
are the same. 

Next we address the physical origin of the coupling term — A0'0. We ffist note that 
this quadratic term is invariant under the exchange of ^ ■0. When A > 0, this term 
can be obtained from a (0 — ^z^)^ term Q, [l^, which represents a local energy penalty 



when the upper and lower monolayers have different compositions. In the case of mixed 
lipid bilayers, such a coupling may result from the conformational confinement of the lipid 
chains, and hence would have entropic origin jl^. By estimating the degree of the lipid 
chain interdigitation, the magnitude of the coupling parameter A was recently estimated by 
May [23:]. In general, the coupling constant A can also be negative depending on the specific 



coupling mechanism 



23| . However, it will be explained later that the phase diagram for 



A < can easily be obtained from the A > one. Hence, it is sufficient to consider only the 
A > case without loss of generality. Although the microscopic origin of the coupling may 
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FIG. 2: (a) Mean-field phase diagram of a single monolayer with a modulated structure in the vicinity 
of the critical temperature, computed using a model as in Eqs. ([5|) - (f7|) . (po is the average composition 
and r is the reduced temperature. The three phases are: striped (S), hexagonal (H), and disordered (D). 
These phases are separated by first-order transition lines, while for clarity we omit showing coexistence 
regions. The filled circle located at (^O)''") = (0,1) indicates the critical point. Note the shift of the 
critical temperature from zero to unity when the modulated phases are considered, (b) Mean-field 
phase diagram of decoupled (A = 0) modulated monolayers at r = 0.8. (po and ipo are the average 
compositions in the two leaflets. The notations of the different phases are described in the text (see 
Sec, nil A]) . All the phases are separated by first-order transition lines. 

differ between systems, we will regard A as a phenomenological parameter and investigate 
its role on the structure, phase behavior and dynamics of coupled modulated bilayers. 

The phase behavior for uncoupled case, A = 0, can be obtained from the analysis of 
-^u[0] 111 is only briefly reviewed here (see also Fig. [2]). For a 2D system, the mean- 
field phase diagram can be constructed by comparing the free energies of striped (S) and 
hexagonal (H) phases. In terms of the order parameter, the stripe phase is described by 



tiQ + 2(f)q cos(g*x) 



(2) 



where 0o = (0) is the spatially averaged composition (imposed by the chemical potential 
fi^), and (pq is the amplitude of the g*-mode in the x-direction. Similarly, the composition 



of the hexagonal phase is given by a superposition of three 2D modes of equal magnitude, 

|qi| = <f 

0H(r) = 00 + ^ cos(qi ■ r), (3) 

V <3 i=l 

where 

qi = q*x, 

q2 = Y (y-x + V?>y 
qs = Y (-X - V^y) , (4) 

and qi = 0. In the above, only the most unstable wavenumber q* is used within the 

single-mode approximation. This can be justified for the weak segregation region close to 
the critical point [sl- 

Averaging over one spatial period, we obtain the free energy densities of the striped, 
hexagonal, and disordered phases, respectively 



/s(0o, 0,) = l<Pl + l<Pt + ir-l + 30^)02 + 1^4^ (5) 



/h(0o, 0.) = 2 00 + + 30^)0? + ^000? + (6) 

/d(0o) = ^0^ + \<l>t (7) 

In Fig. [2](a), we reproduce the original phase diagram of Ref. [3, [ill- The striped, hexag- 
onal, and disordered phases are separated by first-order phase-transition lines. Regions of 
two-phase coexistence do exist but are omitted from the figure for clarity sake [24|. Thus, 
the transition lines indicate the locus of points at which the free energies of two different 
phases cross each other, and are not the proper phase boundaries (binodals). The critical 
point (filled circle) is located at (0o,t) = (0, 1). 



III. TWO COUPLED LEAFLETS WITH THE SAME q* 

Having introduced the free energy and explained the phase behavior of the uncoupled 
case, we shall now explore the equilibrium and non-equilibrium properties of two coupled 
modulated monolayers, A 7^ 0. 
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A. Free Energy Densities 



First we consider the case when D = C = 1 so that the preferred wavenumbers are the 
same for both monolayers, = = q* = ^/ V^- The mean-field phase diagram is calculated 
within the single-mode approximation. Various combinations of 2D modulated structures 
appearing in the two monolayers are possible. The first example is the striped-striped (SS) 
phase, in which both monolayers exhibit the striped phase. This can be expressed as 

0s(r) = 00 + 205 cos(g*a;), (8) 

V's(r) =i>o + 2V'g cos(g*x), (9) 

where 0o = (0) and ipo = {ip) are the average compositions, 0, and ipg are the respective 
amplitudes. These composition profiles are substituted into the free energy of Eq. ([T]). 
Averaging over one spatial period, we obtain the free energy density of the SS phase: 

fss = /s(0o,05) + /slV-o,^,) - A(0o^o + 20,^,), (10) 

where /s is defined in Eq. ([5]). We then minimize fss with respect to both 0^ and ^pq for 
given 00, ipo, t and A. When either 0^ or ipq vanishes, the corresponding monolayer is in its 
disordered phase and the mixed bilayer state will be called the striped-disordered (SD) or 
the disordered-striped (DS) phase. Note that we use the convention that the first index is 
of the 0-leafiet and the second of the -^-one. When both 0^ and ■0^ are zero, the free energy 
density of the disordered-disordered (DD) phase is given by 

/dd = /d(0o) + /d(V^o) - A0oV^o, (11) 



where /d is defined in Eq. (jTl). This free energy /dd was analyzed in Ref. [15|] in order 
to investigate the macro-phase separation of a bilayer membrane with coupled monolayers. 
It was shown that the bilayer can exist in four different phases, and can also exhibit a 
three-phase coexistence. 

Similar to the stripe case, the order parameters of the hexagonal-hexagonal (HH) phase 
can be represented as 

0H(r)=0o + ^X^cos(q,-r), (12) 

^ i=l 

^H(r) = i^o + ■ ^'>- (^^) 

* i=l 



Where the basis of the three qj was defined in Eq. (jlj). By repeating the same procedure as 
for the SS phase, the free energy density of the HH phase is obtained as 



/hh = /h(0O, 0g) + /h(^0, ^5) - A(0O^O + 20^^^^), (14) 

where /h is defined in Eq. ([6]). When either (j)q or ipg vanishes, one of the monolayers is in 
the disordered phase and the bilayer will be called the hexagonal-disordered (HD) phase or 
the disordered-hexagonal (DH) phase. 

When the normal hexagonal phase in one leafiet is coupled to the inverted hexagonal 
phase in the other leafiet, it is energetically favorable to have a particular phase shift of 
27r/3 between the two hexagonal structures. The order parameters which represent such a 
different type of hexagonal-hexagonal (HH*) phase can be written as 

^ i=l 

V^H*(r) = ^o + ^X^ cos(qi-r + y). (16) 

^ i=l 

The free energy density of the HH* phase is then obtained as 

/hh* = /h(0O, 0g) + fnii'O, i'q) - HMo - (l)qi'q)- (17) 

Another combination which should be considered in the present model is the asymmetric 
case where one monolayer exhibits the striped phase and the other the hexagonal phase. 
This striped-hexagonal (SH) phase is expressed as 

0s(r) = 0o + 20gcos(g*x), (18) 



2^ 

VI 

The free energy density of this SH phase is calculated to be 



^H(r) = V^o + -^X]^°^(q^-^)- (19) 

* i=l 



/sH = /s(0o, 0,) + /h(V'o, i^q) - A (^Mo + • (20) 

The phase in which (ps and ip-^ in Eqs. f|T8l) and f|T9l) are interchanged with (/)h and ips is 
called the hexagonal-striped (HS) phase, and its free energy is obtained from the SH phase 
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by noting the cj) ^ ip symmetry. In addition to these phases, we have also taken into account 
the square-square (QQ) phase expressed by 

0Q(r) = 00 + ^ [cos(g*a;) + cos(g*y)] , (21) 

7/>Q(r) = ^0 + ^ [cos(g*x) + cos(g*y)] . (22) 
Then its free energy density is given by 

/qq = /q(</'o, <i>q) + /q(V'o, i^q) - A {Mg + 2(/>,Vg) , (23) 

where 



'0, <Pq} = + \K + ir-l + 3<Pl)<Pl + (24) 

However, we will show below that this QQ phase cannot be more stable than the other 
phases. 



B. Bilayer Phase Diagrams 

Minimizing Eqs. ([10]), ([14]), ([17!), ([20]) and ([2^ with respect to both 0g and tpg, we obtain 
the phase diagram for the coupled bilayer. As a reference, we first show in Fig. [2|(b) the 
phase diagram in the decoupled case (A = 0) for r = 0.8. This can easily be obtained from 
Fig. [2|(a) by combining its two cross-sections (one for 0o and one for i/jq) at r = 0.8. Figure [3] 
gives the phase diagram for a coupled bilayer when (a) A = 0.02 and (b) A = 0.2, while 
the temperature is fixed to r = 0.8 as before. On the {(j)Q,ipQ)-pla:ne, we have identified 
the phase which has the lowest energy, whereas possible phase coexistence regions between 
different phases have been ignored. All the boundary lines indicate first-order transitions. 
Since the free energy Eq. ([1]) is invariant under the exchange of (p ^ ip, the phase diagrams 
are symmetric about the diagonal line 0o = "^o as the upper and lower leaflets have been 
chosen arbitrarily. These phase diagrams are also symmetric under the rotation of 180 
degrees around the origin because Eq. ([1]) is invariant (except the linear terms) under the 
simultaneous transformations of ^ —(f) and ip —>■ This is reasonable as the labels 
of "A" or "B" for the two lipids have been assigned arbitrarily. As a consequence, the 
phase diagrams are also symmetric about the diagonal line 0o = ~^o- The symmetries with 
respect to both 0o = and ■i/'o = in Fig. Mjo) for A = are now broken because of the 
coupling between the two leaflets. 

10 




FIG. 3: Mean-field phase diagram of coupled modulated bilayers for r = 0.8. (/)o and ipo are the 
average compositions in the two leaflets. The coupling parameter is chosen to be (a) A = 0.02 and (b) 
A = 0.2. The notations of the different phases are described in the text (see Sec. Nil A[) . All the phases 
are separated by first-order transition lines. The phase diagram is symmetric with respect to the two 
principal diagonals 0o = V'o and (/)o = — V'o. as described in the text. 

When the couphng parameter is small (A = 0.02), the global topology of the phase 
diagram resembles that of the uncoupled case presented in Fig. Mjo)- Close to the origin, 
00 = ^0 = 0, there is a region of SS phase surrounded by eight other phases: two SH, two 
HS, two HH, and two HH* phases. The HH phase appearing in the region of 0o < and 
V'o < is the combination of the two inverted hexagonal structures on each monolayer. One 
sees that the HH* phase appears in the regions of 0o^o < 0, where the hexagonal and the 
inverted hexagonal structures are coupled to each other. 

A remarkable feature of this phase diagram is the existence of the SS and HH phases 
in the regions where either |0o| or are large. These outer SS and HH phases extend 
up to the maximum or the minimum values of the compositions. These regions of the SS 
and HH phases with A > roughly correspond to those of the SD (DS) and HD (DH) 
phases, respectively, in Fig. EJ^b) with A = 0. Hence the modulated structure in one of the 
monolayers induces the same modulated phase in the other monolayer due to the coupling 
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term. Notice that the SD (DS) phase and HD (DH) phase do not exist in Fig. [3](a). We 
further remark that the extent of the four DD phase regions is almost unaffected by the 
couphng. Even when the temperature is lowered by decreasing r, only the phases located 
close to the origin (0o = V'o = 0) would expand, and the global topology does not change 
substantially. 

For a larger value of the coupling parameter (A = 0.2), the five regions of the SS phase 
merge together forming one single continuous SS region. The four HH regions and still 
distinct and separate the SS region from four DD phase regions. Note that in Fig. Mjo), all 
phases have a symmetric combination of phase modulation such as SS or HH. The asymmet- 
ric combination such as the SH phase does not appear, because the large coupling parameter 
strongly prefers symmetric phases of equal modulations in the two monolayers, although the 
(pq and ipq amplitudes of the two modulated monolayers are not the same in the stripe SS 
phase (or the hexagonal HH phase). As the value of A is increased from 0.02 to 0.2, first 
the SH phase disappears, followed by the disappearance of the HH* phase. When the value 
of A is further increased, the regions of the SS and HH phases expand on the expense of 
the DD phase regions. This means that the coupling between the monolayers causes more 
structural order in the bilayer. Finally we remark that the QQ phase was never found to be 
more stable than any of the other phases considered above. 

Although we have so far assumed that A is positive, the phase diagrams for A < can 
be easily obtained from those for A > by rotating them by 90 degrees around the origin. 
This is because the free energy Eq. ([T]) is invariant under the simultaneous transformations 
of either (p — >■ —0 and A —A, or ■?/' — > —ijj and A — >■ —A. 



C. Modulated Bilayer Dynamics 

In order to check the validity of the obtained phase diagram and to investigate the 
dynamics of coupled modulated bilayers, we consider now the time evolution of the coupled 
equations for and ip- 

Here we have assumed that both and ip are conserved order parameters in each of the 



monolayer (model B in the Hohenberg-Halperin classification 251]). For simplicity, the kinetic 



coefficients and are taken to be unity, and both the hydrodynamic effect and thermal 
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FIG. 4: Equilibrium patterns of coupled modulated monolayers with r = 0.8. In addition setting 
D = C = I implies the same q* in both monolayers. The patterns of (p, ip, (j) + ip, and (p — ip at 
t = 5,000 are presented. The other parameters are chosen to be (a) (pQ = 0.2, ipQ = —0.2, A = 0.02, 
and (b) 0o = 0.5, Vo = 0, A = 0.2. 

fluctuations are neglected. We solve the above equations numerically in 2D using the periodic 
boundary condition. Each simulation starts from a disordered state with a small random 
noise around the average compositions (po and ipQ. In Fig. HI we show typical equilibrium 
patterns of 0, (p+ip and (p—ip for two choices of parameters. The (p+ip pattern is presented 
here because this quantity can be directly observed in the experiment on Montal-Mueller 
bilayers using fluorescence microscopy [12]. The quantity (p — ip measures the concentration 
contrast between the cp and ip leaflets. Time is measured in discrete time steps, and t = 5, 000 
corresponds to a well equilibrated system. In all the simulations below, the temperature is 
flxed to be r = 0.8 corresponding to the weak segregation regime. Notice that all the 
patterns in Fig. H] are presented with the same gray scale. 

Figure Hl^a) illustrates the coupling between a hexagonal phase with 0o = 0.2 and an 
inverted hexagonal phase with ipQ = —0.2 in the weak coupling regime (A = 0.02). Being 
consistent with the phase diagram of Fig. [3]^a), this parameter choice yields the HH* phase 
as seen from the pattern of (p + ip where the two hexagonal structures are superimposed. We 
note that the difference in the order parameter (p — ip also exhibits a hexagonal structure. 



13 



Figure Hl^b) shows the equihbrium patterns for 0o = 0.5 and '?/'o = in the strong 
couphng regime (A = 0.2). If there were no couphng, the 0-monolayer would not exhibit 
any modulation (as it is in its own disordered phase), whereas the t/^- monolayer is in the 
striped phase. We clearly see that, due to the coupling effect, the stripe structure is induced 
in the pattern of 0. This corresponds to the SS phase shown in Fig. [3](b). The periodicities 
of the two striped structures are the same, although their amplitudes differ. Notice that the 
modulation phase of + ^/^ is shifted by vr relatively to that of (p — tp. Since the patterns in 
Fig. m would correspond to the equilibrium configurations, they can be compared with the 
phase diagrams in Fig. [3l We conclude that these simulation results indeed reproduce the 
predicted equilibrium modulated structures. 

IV. COUPLED MONOLAYERS WITH TWO DIFFERENT q* 

We consider next the more general case in which the preferred wavelengths of modulation 
in the two uncoupled leaflets are different, ^ q^. The free energy densities cannot be 
obtained analytically as was done in Sec. IIII Al because there is not a single periodicity on 
which one can average 0(r) and Due to such a difficulty in the analytical treatment, we 
present below the results of numerical simulations, relying on Eq. (!25l) for the time evolution 
of the two coupled order parameters. 

In Fig. [5], we show the patterns for r = 0.8 and 00 = = 0, when both monolayers 
exhibit the striped phase without the coupling. As a reference, we show in Fig. M^a) the case 
when D = C = 1 and A = 0.02 corresponding to the SS phase in Fig. [3](a). The patterns 
of and ip match each other as the composition difference — "0 vanishes throughout the 
system. In Fig. [5](b), the parameters are chosen to be = 0.1296, C = 0.36 and A = 0.02. 
The preferred wavenumbers of the two monolayers are different: q^ = 1.67g^ for uncoupled 
leaflets. The above set of parameters, especially D and C, are chosen in such a way that the 
amplitudes of the two stripes are nearly equal. As long as the coupling parameter is small, 
the two stripes of different periodicities are formed rather independently. The superposition 
of the two striped structures produces an interference pattern resulting in a new modulation 
as seen from the pattern of + "0. The striped modulations of and ip are almost parallel 
or perpendicular to each other. 

When A is made larger, up to A = 0.2 as in Fig. [5](c), the ip field exhibits a complex 
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FIG. 5: Patterns of coupled modulated monolayers with r = 0.8. The patterns of cf), ip, cf) + ip, 
and 4> — are presented for simulation time, t = 5,000, and the average compositions are set to be 
= -^o = 0. The other parameters are chosen to be (a) I) = C = 1, A = 0.02, (b) D = 0.1296, 
C = 0.36, A = 0.02, (c) D = 0.1296, C = 0.36, A = 0.2, and (d) D = 0.1296, C = 0.36, A = 0.4. In 
all cases but (a), D ^ C and the two periodicities are nonequal, q*^ / g^. 

pattern in which two different length scales coexist (reflecting q*^ and g^), whereas the 
pattern of (j) is characterized by a single mode (reflecting q*^). The patters of (f) and ij) 
almost match each other when A = 0.4 as seen in Fig. [5](d). In this case, the modulation 
with a longer wavelength (g^) dominates both monolayers. Figure [5t^b), (c), (d) provide a 
typical sequence of morphological changes, i.e., interference pattern — two-mode pattern — > 
single-mode pattern, as the coupling constant A is increased. 
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FIG. 6: Time evolutions of A defined by Eq. ([26]) . The solid, dashed, dotted, and dashed-dotted lines 
correspond to the parameters of Fig. [5]Ja), (b), (c) and (d), respectively. 

To further analyze the temporal correlations of the two order parameters, and ip, we 
have plotted in Fig. [6] the time evolution of the quantity 



where L = 128 is the linear system size in the simulations. The solid, dashed, dotted, 
and dashed-dotted lines correspond to the time evolutions of A(t) in Fig. M^a), (b), (c), 
and (d) respectively. The solid line (a) first increases and then approaches to zero since 
the patterns of the two stripes coincide in the late stage. The dashed line (b) increases 
in two separate stages. At first because the growth rate of the modulation having smaller 
wavelength (corresponding to ip) is faster than that having larger periodicity (corresponding 
to (p) as is also revealed from the linear stability analysis of Eq. fl2^ which will be published 
elsewhere. At late temporal stages, the value of A remains large for small coupling parameter 
(A = 0.02). When the coupling becomes even stronger (A = 0.2) as for the dotted line (c), 
the value of A is suppressed compared to the dashed line (b), because and ip tend to have 
more overlap for larger A. The same applies for the dashed-dotted line (d) with A = 0.4 as 
compared to the dotted line (c). 




(26) 
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V. DISCUSSION 



We propose a minimal model describing the coupling phenomena between two modu- 
lated bilayers. Considering 2D case, we obtain the mean-field phase diagram when the two 
coupled and spatially modulated monolayers have the same preferred periodicity. Various 
combinations of modulated phases can exist such as the SS, HH, HH* and SH (HS) phases 
as described in Sec. IIIII We have seen that modulations in one of the monolayers induces 
a similar modulations in the other. The region of the induced modulated phase expands as 
the coupling parameter becomes larger. 

When the two monolayers have different inherent wavelengths in the decoupled case, we 
have conducted numerical simulations to investigate the morphologies and dynamics of the 
coupled system. We obtain several complex patterns arising from the frustration induced 
by the two incommensurate structures. As the coupling constant A is made larger, the 
two different modes start to interfere with each other and eventually coincide. The time 
evolution of the striped structures can take place in two steps reflecting the different growth 
rate of the two modulations. 

It is instructive to rewrite the free energy Eq. ([T]) in terms of the sum and the difference 
of and ip with ri± = (p When D = C = 1, we obtain 



F[n+,v-] = j 



(27) 



where /i± = (yU<^ ± /U^)/2. Hence the coupling term between rj^ and //_ in the free energy 
takes the form of an r/^r/^ term with a numerical positive coefficient. The original coupling 
parameter A enters in the coefficients of rj'^ and 77^ terms (but not in the 77^77^ coupling term). 
It shifts the respective transition temperatures of rj^ and ?7_ in opposite directions. When 
the gradient terms are absent and = 0, a similar model was considered by MacKintosh 
and Safran who studied transitions between lamellar and vesicle phases in two-component 



fluid bilayers 



26|. 



The proposed free energy Eq. ([T]) has some analogies to the previous model for the rippled 



phase in hpid bilayers 



271]. It was argued that the coupling between the membrane curvature 



and the asymmetry in the area per molecule between the two monolayers would induce a 



structural modulation of a bilayer. By considering a similar mechanism, Kumar et al. 
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28| 



investigated various modulated phases in two- component bilayer membranes. They claimed 
that the phase behavior of two-component bilayers resembles that of three-component mono- 
layer. This is because the three different local combinations of upper/lower composition in 
bilayers (A/B, B/A, and A/ A for excess of A), would correspond to three different types 
of molecules for the monolayer. One of the new aspects in our model is that the preferred 
wavelengths of the two monolayers can, in general, be different from one another leading to 
a frustrated bilayer state. 

We also point out that there are some similarities between coupled modulated structures 
and the problem of atoms adsorbed on a periodic solid substrate. The latter topic has 
been extensively studied within the Frenkel-Kontorova (FK) model which provides a simple 
description of the commensurate- incommensurate transition [29]. Our model and the FK 
model are analogous in the sense that there are two natural length scales whose ratio changes 
as a function of other model parameters. In the FK model, however, these length scales are 
quenched, whereas in our model they are annealed. 

Another related experimental system can be seen for surface-induced ordering in thin film 
of diblock copolymers jsO]. When the surface is periodically patterned, a tilt of the lamellae 
is induced in order to match the surface periodicity. The situation becomes more complex 
if a copolymer melt is confined between two surfaces. An interesting case arises when the 
spacing between the two surfaces is incommensurate with the lamellar periodicity [s^. 

For systems out of equilibrium, spatial resonances and superposition patterns combining 
stripes and/or hexagons were investigated in a reaction-diffusion model with interacting 
Turing modes of different wavelengths [sil- These models were successful in reproducing 
hexagonal superlattice patterns which are known as "black-eyes" . Although the mechanism 
of pattern formation is different than in our model, we observe similar superposition patterns 
as reported in Ref . 31| , such as hexagons on stripes or hexagons on hexagons (not shown in 
this paper). 

A more detailed study of the present model and several interesting extensions will be pub- 
lished elsewhere. One possible extension is to consider vector order parameters describing. 



for example, the molecular tilt for coupled bilayers 



331] . When the two order parameters 



are vectors, the nature of the transitions between different phases can be different, and even 



the square phase may exist in thermodynamic equilibrium 



3^. 
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